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The present manuscript considers the application of the method of the near-edge X-ray absorption 
spectroscopy (NEXAFS) for the investigation of the graphene-based systems (from free-standing 
graphene to the metal-intercalation-like systems). The NEXAFS spectra for the selected systems 
are calculated in the framework of the approach, which includes the effects of the dynamic core-hole 
screening. The presented spectral changes from system to system are analysed with the help of 
the corresponding band-structure calculations. The obtained results are compared with available 
experimental data demonstrating the excellent agreement between theory and experiment. The 
direct correlation between the strength of the graphene interaction with the metallic substrate and 
the spectral distributions (shape and intensities of vr* and a* features in the C K NEXAFS spectra) 
is found that can be taken as a fingerprint for the description of interaction at the graphene/metal 
interface. 



I. INTRODUCTION 

The method of x-ray absorption spectroscopy (XAS) is one of the widely used spectroscopic tools for electronic 
structure investigations in condensed matter physics [1] . It is based on the absorption of the X-ray photons when core 
shell electron is transferred into the unoccupied valence band states (conduction band) above the Fermi level [Ep] 
(Fig.[T]). Therefore one or more jumps (absorption edges) are usually observed in the absorption spectrum. Moreover 
the energy position of any edge is element specific since it coincides with the energy of the corresponding atomic core 
level. Besides, X-ray transitions are controlled by the dipolar selection rules and for well-defined atomic symmetry 
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of core hole the final state angular momenta are selected in XAS. Strong spatial localisation of the initial core shell 
state leads also to the site specific behaviom- of XAS. 

Generally, appearance of the fine structure of XAS of solids can be described in the framework of photoelectron 
scattering model. Near the edge photoelectron energy is low and photoelectron de Broglie wave undergoes multiple 
scattering on neighbouring atoms. This leads to the origin of the reach and intense fine structure, or near edge X-ray 
absorption fine structure, NEXAFS. On the other hand, in the large photoelectron energies region (far from the XAS 
edge) scattering of the photoelectron waves gives rise to rather simple and weak variation of absorption coefficient, 
or extended X-ray absorption fine structure, EXAFS, well described by single scattering theory. It is considered 
that the range of NEXAFS covers up to 50 eV above the absorption edge and the rest high photon energy part of 
XAS extended up to 1000 eV above the threshold belongs to EXAFS. The latter method is very effective tool for 
the short range atomic ordering determination in multiatomic systems including noncrystalline solids. In NEXAFS 
spectroscopy, the topic of the present paper, absorption spectra can be interpreted on the basis of the corresponding 
molecular orbital or band-structure calculations above Ep. This way is equivalent to the interpretation of XAS data 
in the frame of photoelectron waves multiple scattering model, transparent and demonstrative to date. 

To measure spectral distribution of the X-ray absorption coefficient in massive (bulk) solids yield spectroscopy 
is usually used - the dependence on photon energy of the total photocurrent (so called total electron yield, TEY). 
Since external photo-effect bunch, which accompanied X-ray absorption event, consists of primary fast photo- and 
Auger-electrons together with the cloud of slow secondary electrons, the mean drift way of the whole electron bunch is 
govern by the relatively large mean free path of slow secondary electrons [2, . In order to increase surface sensitivity of 
measurements (for the better representation of the surface-related electronic structure) the negative repulsive potential 
can be applied to the grid-electrode placed in front of the detector to suppress slow kinetic energy component in the 
electron bunch distribution. This experimental mode gives so-called partial electron yield (FEY). In the case of FEY 
NEXAFS with the linearly polarized light, the method can be used for the determination of the orbitals' spatial 
orientation. Generally, one can observe the so-called search-light-like effect [S], which can be used for probing the 
quadrupole moment of the local charge around the absorbing atom. In such XAS experiment, the absorption intensity 
associated with a specific molecular orbital final state has a maximum if the electric field vector is aligned parallel 
to the direction of maximum charge or hole density, i. e. along a molecular orbital, and the intensity vanishes if 
the electric field vector is perpendicular to the orbital axis. A detailed description of the angular dependence of 
NEXAFS intensities can be found elsewhere [31 H] . In case of magnetic samples and circularly polarized light the 
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absorption coefficient depends on the relative orientation of photon spin and magnetization direction of the sample. 
The quantitative analysis of the absorption spectra obtained on the magnetic sample with the circularly polarized 
X-rays can be performed with the help of the so-called magneto-optical sum rules for evaluation of spin- and orbital- 
magnetic moments [3 15] . 

Compared to the relative simplicity of the NEXAFS experiment, the theoretical description of X-ray absorption 
process in solid is not that easy. The main problem here is that the excited system with the core hole has to be 
adequately described, which is always a difficult task for the standard density functional theory (DFT) calculations. 
Particularly it is related to graphite and graphene-based systems. 

The NEXAFS spectrum of graphite or graphene can be separated on two energy regions (see Fig. [ij: (i) Is tt* 
contribution originating from transition of Is core electron on Cpz orbitals and (ii) Is — > a* one originating from 
transition of Is core electron on C sp^ hybrid orbitals (THI]. As was recently shown jTO', although the tt* resonance 
was well reproduced in the framework of the core-hole final state approximation, the shape as well as the intensity of 
the (J* resonance were not reproduced. Only the application of the Mahan-Nozieres-De Dominicis (MND) theory of 
dynamical core-hole screening gives an adequate description of the NEXAFS spectrum of graphite [inillT]. It has to 
be noted here that in the double-peak structure of the a* resonance, the first peak (al), which reflects the excitonic 
final state can be reproduced in the final state and MND approaches, and second (j'2 peak is described as due to the 
quickly completely screened initial-state density of states. 

The present manuscript is devoted to the spectroscopic studies of the graphene-metal systems. The properties of 
these objects are the subject of many experimental and theoretical studies which were recently reviewed in several 
works [T^15j . As discussed the graphene/metal interface is an important object of the graphene-related fields of 
research. Beyond the practical utilisation of these objects (ideal ways to prepare the huge layers of graphene, graphene- 
based spin filters and gas sensors, etc.) [TBH5D], there is also a fundamental interest to the nature of the interaction 
between graphene and metal [TH [TS] . As the NEXAFS spectroscopy can be used for the investigation of the empty 
electronic states of material in the vicinity of Ep, it can be utilised as a complimentary to photoemission spectroscopy 
method, which can shed more light on the nature of interaction between graphene and a metal. 

Here we revisited the theoretical description of the absorption spectroscopy for graphene-based systems. We used the 
scheme described in Refs. [TU||TT] and our computer code for the modelling of NEXAFS spectra of the graphene-based 
systems: graphene single- and bilayer on Ni(lll) and graphene/ Al/Ni(l 11) intercalation-like system. Our results, 
starting from the test example (graphene) and for more complex systems, are discussed on the basis of effects in the 
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electronic structure of the studied objects and compared with available experimental results. Very good agreement 
between theory and experiment is found in our analysis. 

II. THEORETICAL AND COMPUTATIONAL DETAILS 
A. Simulation of the X-ray absorption spectra 

The X-ray absorption spectroscopy, being one of the oldest and most efficient spectroscopic tools to investigate 
the electronic structure of materials, involves a very complicated dynamical many-body process, so that an adequate 
theoretical treatment of the experimental data appears to be a very complicated problem. When the X-ray photon is 
absorbed, the electron from the inner shell excites into the unoccupied states of the conduction band where electronic 
states are described by Wannier functions localized on the absorbed atom. The contribution of transitions into the 
states centred at neighbouring atoms is considered to be small enough to be neglected. 

There are three main approaches for the description the X-ray absorption process. Generally accepted are first two 
one-electron ways of NEXAFS interpretation in terms of the initial state or the final state of system during X-ray 
absorption. Less common is the third way in which many-body processes tracked the X-ray transitions are taking 
into account. These ways are different one of other sonly in approximations of core- hole screening: 

• The first approach assumes that the time of the interaction of the X-ray photon with the electronic subsystem 
is much shorter than any typical process in the system (initial-state approximation). 

• The second approach assumes that the photon-electron interaction is adiabatically slow that the system has a 
time for transfer from the ground state into a new one (final-state approximation). 

• The MND theory assumes the creation of the core-hole as fast enough process, but takes into account the real 
dynamics of the response of the electronic subsystem to a core hole creation (MND approximation). 

1. One-electron approximations. 



In the initial-state approach, the first principle calculations of the electronic structure of material are performed 
with VASP code (for details, see Sec. |IIB ). Then, we obtain the projections of the Bloch wave function onto the 
spherical harmonics centred on site a (atomic number Z) and construct the matrix of initial-state Green's functions 
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as following: 



5mn(w) = ^ 



< a, m|k, I >< k, l\a, n > 



(1) 



cj — £"!(•.; — fd + iS 



where m, n are angular momentum for (s, pi, po, and p-i) orbitals, k - wave vector, I - band index, E - energy, /x - 
chemical potential, and 6 - infinitely small value. Then the initial-state NEXAFS spectrum is given by equation: 



where t„ = E < ls|p|n >, p is the dipole moment operator, and E - X-ray electric field. 

In our work we used these approximations as follows. We model the studied system with a large supercell [i.e. 
(6 X 6) in the case of free-standing graphene]. Then we replace one atom by the following atom in the Periodic Table 
- the so-called Z + 1 approximation. The density of states (DOSs) are obtained by means of VASP calculations 
performed for the large supercell and the Bloch wave functions are projected onto spherical harmonics centred on the 
Z + 1-atom site to get an expression for the final state, and on the Z-atom site to get an expression for the initial 
state. 

The final-state NEXAFS spectrum is calculated according to the equation similar to eqn. ([2|, where the Bloch wave 
fmictions are projected onto atom /3 with an atomic number (Z + 1), and the final-state Green's function Gmn(w) is 
written in the same manner as eqn. (nj). 



The MND theory considers the creation of core-hole as an instant process taking into account real dynamics of 
the response of the electronic subsystems. As a result many-electrons effects due to correlations of electrons with 
core-hole are included. Therefore we construct initial- and final-state Green's functions as mentioned above and the 
so-called core-potential: 




(2) 



2. MND approximation. 



mn 




The MND spectrum is calculated according to the formula shown in the lower part of Fig. [2j The formulas for 



the calculation of the two-particles Green's function [ip{E,t)] and core-hole finite lifetime [A(t)] are presented in the 



left-hand side of the scheme in Fig. [2] The full description of the calculation algorithm can be found in Refs. [Till [TT| 



and in the Supplementary material [26] 
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B. DFT calculations 

The DFT calculations were carried out using the projector augmented wave method "ST], a plane wave basis set 
and the generalized gradient approximation as parameterized by Perdew et al. |22j . as implemented in the VASP 
program [23 . The plane wave kinetic energy cutoff was set to 500 eV. The long-range van der Waals interactions 
were accounted for by means of a semiempirical DFT-D2 approach proposed by Grimme |24] . The corresponding 
structures of the graphene-metal based systems are shown in Fig. [3] and they are discussed in details further in the 
text. The supercell used to model the graphene-metal interface is constructed from a slab of 13 layers of nickel atoms 
with a graphene layer (or bilayer) adsorbed at both sides and a vacuum region of approximately 18 A. In the case of 
graphene/Al/Ni(lll) intercalation-like system, the supercell has (2 x 2) lateral periodicity and contains by one Al 
layer (3 atoms each) introduced between graphene and Ni(lll) on both sides of the slab. In optimizing the geometry, 
the positions (z-coordinates) of the carbon atoms as well as those of the top two layers of metal atoms are allowed to 
relax. In the total energy calculations and during the structural relaxations the fc-meshes for sampling the supercell 
Brillouin zone are chosen to be as dense as 24 x 24 and 12 x 12, respectively. 

To obtain an input for the NEXAFS spectra simulations, we employed large supercells of (6 x 6) and (4 x 4) 
periodicity when studying free-standing graphene and graphene/metal systems, respectively. In the latter case the 
number of nickel layers was reduced to 5. In order to have the same position for the Fermi level, both, initial and 
final states DOSs, are taken from the same calculation sets. 

III. RESULTS AND DISCUSSION 

In this part we present the results of the application of the computational procedure discussed in the previous 
section for the studies of a graphene layer on metallic support. We consider several objects in which the strength of 
interaction of graphene with metal is ranged from "strong" [graphene/Ni(lll)] to "weak" [graphene/Al/Ni(lll) and 
graphene-bilayer/Ni( 1 11)]. 

A. Free-standing graphene vs graphene/Ni(lll) 

In order to prove our implementation of the algorithm discussed earlier we performed the calculation of NEXAFS 
spectra for the free standing graphene (Figs. |4] and [s]). The main peaks of the graphene spectrum at 1.9eV and 8.3eV 
above Ep correspond to transitions of the Is core electron onto the tt* and a* states, respectively. Fig. |4] shows the 
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theoretical spectra calculated according to the procedure described in Sec.|ll]for ten different incidence angles a. One 
can note a smooth decreasing (increasing) weights of Is — )• tt* (Is — )• a*) transitions when going from a = 0° to 
a = 90°, that is in accordance with geometry of the absorption process. This example is a nice demonstration of 
the search-light-like effect discussed earlier: for the grazing incident angles the maximum of absorption is observed 
for TT* states and for angles close to the normal to the surface the maximum of absorption is detected for Is — > a* 
transitions. Furthermore, it is clearly seen that both the shape and positions as well as relative intensities of the 
peaks are well reproduced by theory [for comparison between calculated and experimental spectra PHI US] , see Fig. [s] 
(upper panel)]. As was discussed earlier in Refs. pill [TTj the shape as well as the position of Is tt and Is a\ 
transitions are well reproduced by the MND theory described in Fig. [2] The second peak a\ can be reproduced only 
if fully screened initial states density of states is taken into account and the corresponding initial states NEXAFS 
spectrum is summed with MND spectrum in order to reproduce the experimental observations. 

Before investigating modifications in the absorption spectra due to the presence of the Ni(lll) substrate under- 
neath graphene, let us consider briefly the atomic structure of the interface. The surface lattice constant of Ni(lll) 
matches the lattice constant of graphene almost perfectly and therefore (1x1) structure is formed. The structure of 
graphene/Ni(lll) widely accepted in the literature has the carbon atoms arranged in the so-called top-fee configura- 
tion on Ni(lll): one carbon atom from the graphene unit cell is placed above the Ni interface atom and the second one 
occupies the hollow fee site of the Ni(lll) surface (Fig. [s]). Our calculations support this model (see Supplementary 
Tab. SI US]). 

The equilibrium distance between graphene layer and Ni(lll) is relatively small, indicating "strongly" bonded 
system: significant hybridization of the valence band states of graphene and Ni leads to the complete destruction of 
the graphene Dirac cone [121 [?7H^ . [This effect of hybridization is clearly seen when considering electronic band 
structure of the system (see Supplementary Fig. SI [21]) ■ The modifications in the electronic structure are considered 
in details in Ref. p7]]. 

Consequently, the NEXAFS spectra in this case are modified with respect to those of free-standing graphene [Figjs] 
(middle panel)]. Firstly, in the region of the Is — > tt* transition the graphene/Ni(lll) absorption spectra have a 
double-peak structure compared to those of graphite that is explained by the transitions of C Is core electron into 
two unoccupied states (interface states), which are the result of hybridization of the graphene tt and Ni 3d valence 
band states. The first peak at about 0.2 eV corresponds to the interface state which is allocated in the range between 
and 1 cV above Ep around if-point of the Brillouin zone, having slightly different energy depending on the type of 
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C-atom (C^f or 0(^1") and on the spin-orientation. The second peak, near 2 eV, comes mainly from the unoccupied 
interface state at the AI point (see Supplementary Figs. SI and S2 [26J. From this analysis one can clearly see that 
both peaks belong to the hybrid states where both carbon atoms are involved. Therefore the previous analysis of 
NEXAFS spectra for the graphene/Ni(lll) system performed in Refs. [inillSKSn] on the basis of calculations of the 
K-edge electron loss spectra 27J is not fully correct (the double peak structure was assigned to different carbon atoms 
in the unit cell). Secondly, the visible reduction in the energy separation between the tt* and cr* features are observed 
for all graphene/Ni(lll) spectra compared to that in the spectra of free-standing graphene, in accordance with the 
shift of the a* state for the carbon atoms at the interface by ca. 1 eV to lower energies that can be explained by the 
lateral bond softening within the adsorbed graphene monolayer and can be connected with the downward shift of a* 
states for graphene/Ni(lll) compared to free-standing graphene. 

B. Decoupling graphene from the substrate: The graphene/Al/Ni(lll) system 

In the graphene/Al/Ni(lll) system, the fee Al(lll) monolayer is placed in the space between graphene and Ni(lll) 
(Fig. |3]). There is only one possible arrangement of aluminium atoms in the high-symmetry positions as shown in 
the figure: the Al(lll) lattice plane is rotated by 30° with respect to the graphene/N(lll) lattice such that the Al 
atoms occupy all three kind of high-symmetry adsorption sites {top, hollow fee, hollow hep) in the space between the 
graphene layer and the Ni(lll) surface. In this structure there are four different occupation sites for carbon atoms 
in the unit cell with respect to the adsorption sites of the Ni(lll) surface when they are placed above a Ni or Al 
atom. They are labeled in Fig. [s] as C^f , C^'?'^, C^j^, and C;^j^, respectively. A mean distance between the graphene 
and Al layer is 3.312 A, indicating "weak" bonded graphene-metal system [31]. The latter fact can be evaluated when 
considering the electronic band structure of the system under study: all electronic bands of graphene are shifted to 
lower binding energies compared to graphene/Ni(lll). Furthermore, the electronic structure of the graphene layer as 
well as the Dirac cone in the vicinity of Ep are fully restored (there is a small electron doping of graphene leading 
to a shift of the Dirac point below Ep by ca. 0.64 eV). Thus, intercalation of Al decouples the electronic structure of 
graphene from the substrate, preventing thereby hybridization between graphene tt and Ni 3d valence band states. 
In the present system there is no strong hybridization between graphene and Al valence band states conserving the 
Dirac cone. There is only one energy region where such hybridization is visible: around the M point at approximately 
E~Ep^leY (see Ref. [31]). 

The decoupling process in this case is also easily visible when considering the NEXAFS spectra obtained for the 
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system of interest. Intercalation of a thin Al layer underneath a graphene layer on Ni(lll) leads to drastic changes 
in the absorption spectrum compared to the one for graphene/Ni(lll): The shape of the spectrum, positions of main 
spectroscopic features as well as the energy separation between tt* and a* features become similar to those in the 
spectrum of free-standing graphene [Fig. [S] (lower panel)]. Due to the avoiding hybridisation between the valence 
band states of graphene and metal in this case, the individual spectra evaluated for the C-atoms of the four types are 
indistinguishable from each other (see Supplementary Fig. S3 [2B]). Weak hybridisation between graphene tt* and Al 
3p states does not affect the result. Note: agreement between theory and experiment is again excellent. 

C. Bilayer-graphene/Ni(lll) 

The energetically most stable arrangement of the graphene layers above Ni(lll) in this case corresponds to the 
situation when carbon atoms of the first layer are adopting the top-fcc configuration [as it is in the case of the single 
graphene layer deposited on Ni(lll)] and the second graphene layer is displaced horizontally with respect to the first 
graphene layer as it is in bulk graphite in the way, that carbon atoms from the graphene unit cell are occupying the 
both hollow sites with respect to the Ni(lll) surface [Fig. [s] (right-hand side figure)]. The resulting configuration can 
therefore be named as top-fcc/fcc-hcp (see Supplementary Tab. SI [26 ). 

In the studied system, the interaction energy between the top (second) graphene layer and graphene/Ni(lll) is 
91 meV per graphene unit cell that is comparable with the value known for the interlayer interaction in graphite 
(122 ± lOmeV 32 ). Therefore we can expect, that bilayer-graphene/Ni(lll) should display properties of the both 
systems, free-standing graphene and graphene/Ni(lll). When looking at the results of the electronic band structure 
calculations, one can see that, indeed, while valence band states of the first graphene layer and Ni are strongly 
hybridized and Dirac cone is destroyed (see Supplementary Fig. S4), picture observed for the second graphene 
layer is much more similar to that of the free-standing graphene: there is no hybridisation visible and Dirac cone is 
preserved (see Supplementary Fig. S5 [35]). At the same time, a band gap of ca. O.lSeV is opened at the X-point 
at E — Ep = — 0.5eV. There is, however, no significant difference between pictures observed for the empty states of 
the second graphene layer and those of the free-standing graphene. The same is true for the empty states of the first 
graphene layer in this system and graphene in graphene/Ni(lll). This is reflected in the NEXAFS spectra calculated 
for the bilayer: while results obtained for the first graphene layer are identical with those of graphene/Ni(lll), the 
NEXAFS spectra of the second layer are similar to those obtained for the free-standing graphene as well as graphene 
in the graphene/ Al/Ni(l 11) trilayer [Fig. |6] (left panel)]. 
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As discussed earlier the NEXAFS spectra can be collected in two modes, more bulk- sensitive TEY and more 
surface-sensitive FEY modes. Here, the variation of the repulsive potential (which suppress the secondary electron 
tail) applied to the electrode before electron detector allows for the gradual variation of the sensitivity (surface vs 
bulk). Fig. [6] (right panel) shows the NEXAFS spectra of graphene-bilayer/Ni(lll) where the sensitivity is varied 
from TEY mode (lower spectrum) to PEY mode (upper spectra) . The clear change of the shape of the spectra to the 
one characteristic for the free-standing graphene is visible. These simulations allow to discriminate between single 
and bilayer graphene in the investigated system. 

IV. CONCLUSIONS 

Here we presented the theoretical description of the NEXAFS spectra of the graphene layers on metallic substrate. 
The interaction in these systems is varied from "strong" (when the electronic structure of graphene in the vicinity of Ep 
is strongly modified) to "weak" (when the electronic structure of the graphene layer is a free-standing-graphene like) . 
All these variations in interaction and in electronic structure are consequently reflected in the simulated NEXAFS 
spectra which can be used as a fingerprint for the description of interaction at the graphene/metal interface. The 
presented results are compared with the available experimental data and they show the excellent agreement between 
theory and experiment. 
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in-plane 




FIG. 1: The scheme of the X-ray hght absorption process. The hght from the synchrotron Hght source is scanned around the 
energy corresponding to the binding energy of the particular core level and in this case the absorption coefficient is proportional 
to the density of the valence band states above Ef- The orbital selectivity is reached via relative orientation of the linearly 
polarised light and the sample. 
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FIG. 2: Computational algorithm used for calculation of NEXAFS spectra. 
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FIG. 3: Top- and side-views of the crystallographic structures of the studied systems: graphene/Ni(lll), graphene/Al/Ni(lll), 
and graphene-bilayer/Ni(lll) (from left to right). 
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FIG. 4: Calculated, in a framework of MND approach, C K NEXAFS spectra for a free-standing graphene as a function of the 
angle, a, between electric field vector of X-ray and the normal to the surface. The experimental geometry is shown as an inset. 
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FIG. 5: Comparison of experimental and calculated NEXAFS spectra obtained with a — 40° at the C K absorption edge for 
graphite, graphene/Ni(lll), and graphene/Al/Ni(lll). Experimental data are taken from Refs. [25] and [31j . 
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FIG. 6: (Left panel) Calculated C K NEXAFS spectra for a = 40° for the graphene-bilayer/Ni(lll) system. Spectra are 
separated on the contributions from the 1st and 2nd layer, respectively. (Right panel) Sequence of the simulated C K NEXAFS 
spectra for graphene-bilayer/Ni(lll), where the contribution from both layers is increased/decreased from top to bottom in 
order to model the change of surface- vs-bulk sensitivity in experimental data. 
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Content: 

1. Description of the calculation algorithm of the NEXAFS spectra. 

2. Table SI: Results for the graphene/Ni(lll) and graphene-bilayer/Ni(lll). All adsorption energies (AE) are 
given per graphene unit cell. 

3. Figure SI: Calculated majority and minority spin band structures (black) for a slab terminated by 
graphene/Ni(lll) interface {top-fcc arrangement). The bands replotted in red (blue) using the carbon 
{sp"^) character as a weighting factor are superimposed. Contributions from the C*°p and C^'^'^ are separated. 

4. Figure S2: Calculated C K NEXAFS spectra for a = 40° for the graphene/Ni(lll) system. Spectra are 

separated on the contributions from the C*°^ and C^"'^, respectively. 

5. Figure S3: Calculated C K NEXAFS spectra for a = 40° for the graphene/ Al/Ni(l 11) system. Spectra are 
separated on the contributions from the C^f , C^f , C;^;*^, and C/jj^, respectively. 

6. Figure S4: Calculated majority and minority spin band structures (black) for a slab terminated by graphene- 
bilayer/Ni(lll) interface {top-fcc/ fcc-hcp arrangement). The bands replotted in red (blue) using the carbon p^ 
{sp^) character of the first graphene layer as a weighting factor are superimposed. Contributions from the C*"^ 
and C^'^'^ are separated. 

7. Figure S5: Calculated majority and minority spin band structures (black) for a slab terminated by graphene- 
bilayer/Ni(lll) interface [top-fcc/ fcc-hcp arrangement). The bands replotted in red (blue) using the carbon 
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(sp^) character of the second graphene layer as a weighting factor are superimposed. Contributions from the 
Qhcp g^jj^ Qfcc gj.g separated. 

8. Figure S6: Calculated C K NEXAFS spectra for a = 40° for the free-standing graphene as obtained by means 
of different approximations (SK = Slater-Koster method). 

9. Figure S7: Calculated C K NEXAFS spectra for a = 40° for the graphene/Ni(lll) system as obtained by 
means of different approximations (SK = Slater-Koster method). 

10. Figure S8: Calculated C K NEXAFS spectra for a = 40° for the graphene/Al/Ni(lll) system as obtained by 
means of different approximations (SK = Slater-Koster method). 

11. Figure S9: Calculated C K NEXAFS spectra for a = 40° for the graphene-bilayer/Ni(lll) system as obtained 
by means of different approximations (SK = Slater-Koster method). 
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Description of the calculation algorithm of the NEXAFS spectra. 



Let us consider the radiative X-ray transitions of an electron from the core level c = {nc,lc,Tnc} on the band 
states nk localized on the atom A. In the tight-binding approximation, neglecting the "cross" transition contribution 
to the intensity of radiation, the NEXAFS intensity I{E) can be represented as the sum of the products of the 
electric dipole probabilities of transition W/^ from the core level on the atom-like valence band state with partial 

density of states Nn{E) 

/(£)~^W^^7V^(£;), (3) 
ft 

W^^~El^Mmj'> (4) 
where the amplitude A^^^ of transition is given by 

^tm. = i^M'Pc) ■ (5) 

Here and ipc are the atomic wave functions of valence electron and core hole, respectively. In the framework of 
the one-particle approximation, the partial density of states (PDOS) can be expressed via the imaginary part of the 
advanced Green's function projected onto the atomic orbitals G~ (E) 

N^{E) = ^lmG-^{E). (6) 

Completely ignoring the influence of the core hole, the Green's function can be written as 

r<- lT?\-\^ (MlV'nk) (V'»k|'^) 

^ft^i^)-L.E-En{k)-E, + iS' (7) 
with E{ - the Fermi energy, -B„(k) and Vnk - the unperturbed band energy and wave function, respectively. It may 
be considered as an assumption that the lifetime of the core hole is much smaller than the relaxation time of the 
system (an approximation of the final state) . 

Another extreme case takes place when the band wave functions are completely changed because of the core hole 
Coulomb potential (initial state approximation). Then the expression for the Green's function takes the form 

G-AE) = V ^^L^"*-^ ^^"'-'"^ , (8) 
' ^E-Er.{k)-Ef + iS' ^ ' 

where -En(k) and are perturbed band energies and wave functions, respectively. For PDOS in the presence of the 

core hole one gets 

N^{E) = -ImG-^iE). (9) 
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MND theory considers the core hole filHng completely in the dynamical way. The expression for the NEXAFS 
intensity within the MND theory in the form 

oo +00 

lMnY,{E)^ReY,wJ dE' N^,{E') / x (^+^-^') *+^(-*) ^.^(i?', t) , (10) 



where the matrix Lp[E,t) satisfies the system of integral equations 

00 

(^^,(^,t) -(5^,+ ^ V. j dE' X L^,,,{E,E',t)^,,,iE',t). 



(11) 



The kernel of the integral operator L^^{E,E' ,t) is given by 



L,.{E,E',t) = Y. ^^^'^^'''^ |<^^J')%.-(g,t) ^ ^^2) 



where 



I^..{E, t) = 5,,, + J2 y,,' I dE' X N^.^E') , (13) 



-00 



V - ((/^l^»k)(V^„kk)£^n(k) - (/^|^„k)(^«k|i^)^n(k)) . (14) 



In the above expressions, the matrixes N^^, and iV^^ have the form 



N^, = - Im {E) , N^, = - Im G;, (S) . (15) 



Note that the diagonal elements of matrixes A^^i, and N^i, coincide with the PDOS defined above in Eqs. (|6|) and (|9 
in the final and initial state approximations. 



There are two ways to define input data for the MND calculations expressed by Eqs. (10|-(13). The first approach 



is to use the expressions ([7|, ([8|), and (14 1 to introduce the Green functions G^{E),G^{E) and the core potential 
V . It works well if the band wave functions basis set is large enough. If it is not the case, the core hole potential V 
depends on energy. To minimize the error caused by the energy dependence of the potential one can use a second 



approach. The unperturbed Green function G (E) and core hole potential V defined from equations (|8| and (14 1, 
respectively, whereas the perturbed Green function G^ (E) is obtained from the Slater-Koster equation 

g-^{e) = Y^g;Ae)[t-']^,^, (16) 

where 

T,.. = 5^.+Y,V^,,G-,^{E) (17) 
Although the potential still depends on energy, the perturbed Green function now takes it into account. 
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Table SI. Results for the graphene/Ni(lll) and graphene-bilayer/Ni(lll). All adsorption energies (AE) are given per 
graphene unit cell. 



System 



AEi (eV) 



AE2 (eV) 



AEs (eV) 



graphene/Ni(lll) 

top-fcc 

top-hep 
fcc-hcp 



-0.284 

-0.257 
-0.060 



gr aphene-bilayer /Ni(lll) 

top-fcc / top-hep 
top-fcc / fcc-hcp 
top-fcc / top-fcc 



-0.091 
-0.096 
-0.067 



-0.187 
-0.190 
-0.175 



-0.275 
-0.280 
-0.277 



top-hep / top-fcc 
top-hep / fcc-hcp 
top-hep / top-hep 



-0.089 
-0.096 
-0.067 



-0.173 
-0.176 
-0.162 



-0.247 
-0.253 
-0.250 



fcc-hcp / top-hep 
fcc-hcp / top-fee 
fcc-hcp / fcc-hcp 



-0.128 

-0.140 
-0.101 



-0.094 

-0.100 
-0.080 



-0.088 

-0.101 
-0.087 



AEl — -Egr2/Ni(lll) — -E'gr/Ni{lll) " -Egr 

AE2 = i?gr„/Ni(iii) — -ENi(iii) — -Egr X H (ti - number of graphene layers) 

AE3 = -Egr2/Ni(lll) — -E'Ni(lll) — Egr2 




Figure SI. Calculated majority and minority spin band structures (black) for a slab terminated by graphene/Ni(lll) 
interface itop-fcc arrangement). The bands replotted in red (blue) using the carbon (sp^) character as a weighting 
factor are superimposed. Contributions from the C*°p and G^'^'^ are separated. 
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Figure S2. Calculated C K NEXAFS spectra for a = 40° for the graphene/Ni(lll) system. Spectra are separated 
on the contributions from the C*°p and C^'^'^, respectively. 
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Figure S3. Calculated C K NEXAFS spectra for a — 40° for the graphene/Al/Ni(lll) system. Spectra are separated 
on the contributfons from the C^j^, Cj^?, C;^'j^, and C^'?'^, respectively. 
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Figure S4. Calculated majority and minority spin band structures (black) for a slab terminated by graphene- 
bilayer/Ni(lll) interface {top-fcc/fcc-hcp arrangement). The bands replotted in red (blue) using the carbon (sp^) 
character of the first graphene layer as a weighting factor are superimposed. Contributions from the C*°p and C^"^"^ 
are separated. 
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Figure S5. Calculated majority and minority spin band structures (black) for a slab terminated by graphene- 
bilayer/Ni(lll) interface {top-fcc/fcc-hcp arrangement). The bands replotted in red (blue) using the carbon (sp^) 
character of the second graphene layer as a weighting factor are superimposed. Contributions from the C'"^^ and C^"^"^ 
are separated. 
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Figure S6. Calculated C K NEXAFS spectra for a — 40° for the free-standing graphene as obtained by means of 
different approximations (SK = Slater-Koster method). 
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Figure S7. Calculated C K NEXAFS spectra for a = 40° for the graphene/Ni(lll) system as obtained by means of 
different approximations (SK = Slater-Koster method). 
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Figure S8. Calculated C K NEXAFS spectra for a = 40° for the graphene/Al/Ni(lll) system as obtained by means 
of different approximations (SK — Slater- Koster method). 
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Figure S9. Calculated C K NEXAFS spectra for a — 40° for the graphene-bilayer/Ni(lll) system as obtained by 
means of different approximations (SK — Slater-Koster method). 



